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Biomass torrefaction: Modeling of volatile and solid product evolution kinetics 
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Model developed to describe evolution of volatiles during willow torrefaction. 

► Composition of volatiles defined in terms of nine identifiable chemicals. 

► Solid residue composition defined in terms of carbon, hydrogen and oxygen content. 

► Coupled to solid mass loss kinetic rate equations. 
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The aim of this work is the development of a kinetics model for the evolution of the volatile and solid 
product composition during torrefaction conditions between 200 and 300 °C. Coupled to an existing 
two step solid mass loss kinetics mechanism, this model describes the volatile release kinetics in terms 
of a set of identifiable chemical components, permitting the solid product composition to be estimated by 
mass conservation. Results show that most of the volatiles released during the first stage include highly 
oxygenated species such as water, acetic acid, and carbon dioxide, while volatiles released during the sec¬ 
ond step are composed primarily of lactic acid, methanol, and acetic acid. This kinetics model will be used 
in the development of a model to describe reaction energy balance and heat release dynamics. 

© 2012 Elsevier Ltd. All rights reserved. 


1. Introduction 

Torrefaction is a thermochemical pretreatment process which 
improves biomass utilization characteristics such as gravimetric 
heating value, grindability as measured by the grinding energy 
per unit heating value, and hydrophobicity (Acharjee et al., 2011; 
van der Stelt et al., 2011 ). During torrefaction, biomass is treated 
at temperatures between 200 and 300 °C, in an inert atmosphere 
for a period usually between 15 and 60 min. 

The relevant process parameters during torrefaction include the 
temperature, reaction time, heating rate, feedstock, and particle 
size. The chemical characteristics - including ultimate analysis - 
of biomass torrefied under numerous combinations of tempera¬ 
ture, reaction time, and feedstock have been investigated experi¬ 
mentally by Prins (2005), Arias et al. (2008), Bridgeman et al. 
(2008), Almeida et al. (2010) and Medic et al. (2011). Typically, 
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the torrefied product retains 80-95% of the energy and 70-90% 
of the mass of the original raw biomass. The remaining 10-30% 
of the mass is released in the form of gaseous species. Prins 
(2005) and van der Stelt (2011) have analyzed the evolution of vol¬ 
atile species released during torrefaction utilizing gas chromatog¬ 
raphy (GC) and high performance liquid chromatography (HPLC). 

While experimental work on torrefaction remains active, few 
models exist to comprehensively describe the evolution of the vol¬ 
atile products or solid composition and energy balance over a 
range of conditions. To address this need, we have developed a ser¬ 
ies of models to assist in the study of torrefaction according to the 
following steps: 

1. A kinetics model to describe the release of volatiles and the 
resulting change in composition of the solid product. 

2. A thermochemical model to estimate product properties (e.g. 
specific heat and heat of formation) and reaction energy balance 
based on the changes in chemical composition predicted by the 
volatile release model. 

3. Coupling of this thermochemical model to a single particle heat 
and mass transfer model of torrefaction to study process condi¬ 
tions and particle size effects. 
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This paper describes the kinetics model of solid and volatile 
product evolution, which was used as the basis for the subsequent 
models. Model predictions of solid and volatile product character¬ 
istics were compared with published experimental data. Models 
for the thermochemistry and single particle torrefaction will be 
published in follow up articles. 

Models for the evolution of the solid product characteristics 
have been developed based on experimental characterization un¬ 
der specific torrefaction conditions and feedstocks. One kinetics 
model proposed and validated by Prins (2005) describes the solid 
mass loss of willow during torrefaction. Based on corn stover tor- 
refaction experiments, (Medic et al„ 2011) fitted quadratic regres¬ 
sions relating process parameters such as initial moisture content, 
temperature, and residence time to solid product energy and mass 
yield. Based on eucalyptus torrefaction experiments (Almeida 
et al„ 2010) showed that solid mass loss can be used as a quantita¬ 
tive indicator of the extent of torrefaction. Several linear regres¬ 
sions were developed relating the mass loss to solid product 
characteristics such as the energy yield, heating value, and chemi¬ 
cal composition). These empirical regressions do not describe the 
composition of the volatile or solid torrefaction products. Using 
the static dessicator technique, the effect of torrefaction on the 
equilibrium moisture content of raw loblolly pine was represented 
by four empirical parameters (Acharjee et al., 2011). Increased deg¬ 
radation of the hemicellulose component of biomass was found to 
correspond with an increasingly hydrophobic solid product. 

A detailed multi-step chemical kinetic model of pyrolysis has 
been developed by Ranzi et al. (2008). The validity of this long 
mechanism under torrefaction conditions, which is associated with 
slow heating rates and lower temperatures, is still underway. 
Neves et al. (2011) developed an empirical model for biomass 
pyrolysis products based on fitting the trends of experimental data 
between (200 and 1000 °C). However, of the seven modeled vola¬ 
tile constituents, some pyrolysis products (“C x H y ”) are not defined 
as readily identifiable chemicals. Additionally, the chemical com¬ 
position of the resulting char depends only on a single parameter, 
the peak pyrolysis temperature. Therefore, this empirical model is 
unable to represent the effect of residence time on solid product 
composition during torrefaction. 


2. Model 

2.1. Approach 

The present approach built upon the existing mass loss kinetics 
model for torrefaction published in (Prins, 2005). This kinetics 
model enables accurate prediction of the solid mass balance during 
torrefaction under kinetically limited conditions in the absence of 
transport limitations. Within the kinetics mechanism, solid and 
volatile products are lumped into solid (A, B, C) and volatile (1/1, 
V2) pseudocomponents; however, their composition in terms of 
actual chemical species and the enthalpy of reaction for each step 
has not been estimated or modeled previously. 

Separate volatile composition experiments using HPLC also by 
Prins (2005) reported the cumulative yield of nine chemical species 
including acetic acid, water, formic acid, methanol, lactic acid, 
furfural, hydroxyacetone (acetol), carbon dioxide and carbon mon¬ 
oxide which were generated during willow torrefaction under five 
torrefaction conditions. 

The present approach to developing a simplified volatile com¬ 
position model was comprised of three steps: 

(1) Analysis of appropriate solid mass loss kinetics mechanism 
and volatile composition data for torrefaction from the liter¬ 
ature (Section 2.2 and 2.3). 


(2) Application of fitting techniques to published experimental 
data to enable the formulation of a volatile composition 
model describing the volatile pseudo components (1/1 and 
V2) in terms of nine major chemical components (acetic 
acid, water, formic acid, etc.) (Section 2.4). 

(3) The application of mass conservation to enable formulation 
of a solid composition model describing the ultimate analy¬ 
sis of the solid product and its constituent pseudo compo¬ 
nents, B and C (Section 2.5). 

2.2. Solid mass loss kinetics model 

Prins showed that a competitive, two stage, first order mecha¬ 
nism accurately modeled the solid mass loss profiles during kinet¬ 
ically controlled torrefaction conditions of willow: 




C 


where A is the raw biomass (willow), B is a solid intermediate reac¬ 
tion product, C is a residual solid product, and VI and V2 represent 
volatiles. This mechanism was originally proposed by Di Blasi and 
Lanzetta (1997) to describe pure hemicellulose decomposition. 
The four Arrhenius kinetic parameters (ki, k 2 , k vu fc^) were fitted 
to experimentally measured mass loss data. At any moment in time, 
the mass of the solid product is described by the sum of masses of A, 
B, and C, while the total mass of volatiles is provided by the sum of 
VI and V2. The equations describe the evolution of these pseudo 
components: 


dm A ,, , . 

r A = ~^r=-(ki + kvi)*m A 

0) 

dm B . . 

r„ = -gp = (fei) *m A - (k 2 + few) * m B 

(2) 

dm c „ , 
r c = ^ = (k 2 )*m B 

(3) 

r vi = ^p = (fcvi) * m A 

(4) 

r v 2 = -gp = (kvz) * m B 

(5) 

where m x is the mass of pseudo component (X = A, B, C, VI, V2), 
r x = represent the net mass production rate, (the net rate of 
change of mass of the pseudo component X) in (kgXs ’). The 
parameters that Prins’ fitted during willow TGA mass loss experi¬ 
ments are: 

fc 4 = 2.48 x 10 4 exp(-75976/RT) 

(6) 

few =3.23 x 10 7 exp(-114214/RT) 

(7) 

fe 2 = 1.1 x 10 10 exp(-151711/RT) 

(8) 

few = 1.59 x 10 10 exp(-151711/RT) 

(9) 


where the preexponential factor has units of s ', the activation 
energy has units of J mol \ R is the universal gas constant in 
J moU 1 K ', and T is the particle temperature. 

As was noted in Prins (2005) the first reaction phase 
(A -> B + VI) is much faster than the second phase (B-> C+V2). 
Solid mass loss kinetic experiments and analyses by van der Stelt 
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(2011) also concluded that torrefaction can be modeled in two 
stages. The dynamics of this mechanism can be explained by the 
pyrolysis behavior of pure lignocellulose components (hemicellu- 
lose, cellulose, and lignin) summarized next: 

Hemicellulose, cellulose, and lignin comprise 17-20%, 53-60%, 
and 25%, respectively of a hardwood such as willow (Bergman 
et al„ 2005; Bridgeman et al„ 2008). During TGA mass loss exper¬ 
iments of pure components, hemicellulose has been found to be 
the most reactive component, decomposing between 200 and 
300 °C while cellulose degrades between 275 and 350 °C (Biagini 
et al„ 2006). Lignin is the least reactive and decomposes over the 
range of 200 and 600 °C. 

Simulation of the solid mass loss (using the willow kinetic 
parameters) indicate that the first stage concludes within 15- 
60 min, and results in 16-30% solid mass loss for temperatures be¬ 
tween 250 and 300 °C. The second stage takes up to several hours 
to reach completion and results in an additional 42-48% mass loss; 
thus, the mass loss occurring during the faster first stage of torre¬ 
faction is primarily attributable to the decomposition of hemicellu¬ 
lose (with an increasing contribution from cellulose decomposition 
at higher temperatures). The mass loss during the slower second 
stage is primarily due to cellulose decomposition, with minor lig¬ 
nin decomposition and charring of the remaining hemicellulose. 

Recent experimental analysis also supports this description. The 
torrefaction of loblolly pine at three levels of severity (16.5%, 24.1%, 
and 46.5% mass loss) demonstrated that the hemicellulose experi¬ 
enced almost complete degradation under the most severe condi¬ 
tion (Meng et al„ 2012). The cellulose experienced only partial 
decomposition while the lignin fraction remained mostly intact. 

Though the mechanistic description of torrefaction in terms of 
the pyrolysis of its pure lignocellulose components was originally 
propounded by several authors Bergman et al. (2005) and Prins 
(2005) it has been supported by further TGA experiments (Chen 
and Kuo, 2011; van der Stelt, 2011 ). 

2.3. Volatile composition experimental data 

In a separate analysis from his mass loss kinetics parameter fit¬ 
ting, Prins reported the average composition and mass of volatiles 
produced during five different experimental conditions during wil¬ 
low torrefaction (see Table 1) with solid mass yields between 
66.7% and 92.6%. The experimental conditions as well as the cumu¬ 
lative yields of nine different chemical species collected are 


reported in Table 1. Batch experiments of 10 grams of biomass 
with particle diameters between 0.7 and 2.0 mm accompanied by 
slow heating rates of 10°Cmin 1 ensured that the torrefaction 
was kinetically controlled (Prins, 2005). The condensable volatile 
products were liquefied in a cold trap at 5 °C while non-condens- 
able products were collected in a gas bag. All products were 
weighed to produce mass balance and the volatiles were analyzed 
using gas chromatography and high performance liquid 
chromatography. 

Torrefaction reactions initiate at 200 °C. T final is the maximum 
temperature the reactor reaches after the heating period. The hea¬ 
tup period, which occurs between when the temperature is ramp¬ 
ing from 200 °C to Tfi n ai, is determined by the heating rate. The total 
time is the summation of the heatup period and the isothermal 
period (the time the sample is held at T fina i). The total time does 
not include the time to heat the sample from ambient temperature 
to 200 °C, nor does it include the time when the particle is being 
cooled. The solid and volatile mass yields are defined in terms of 
the percent weight of the initial solid mass, which is taken when 
the temperature is 200 °C. 

It is clear from the experimental results in Table 1 that the aver¬ 
age composition of the total volatile yields changes with torrefac¬ 
tion temperature, residence time, and mass yield. For example in 
experiment 1, methanol comprises 2.4% of the total volatile yield 
while in experiment 5 methanol comprises 11% of the total volatile 
yield. 

In order to verify the assumption that the experimentally mea¬ 
sured volatiles did not undergo significant intra particle or extra 
particle secondary decomposition, a characteristic time associated 
with the cracking of condensible volatile products was calculated 
from experimentally measured Arrhenius kinetic data for beech 
wood (Rath et al., 2002): 

W = 1 /Ae- E,RT (10) 

where A = 10 514 (s ')and E = 93.37 (kj mor 1 ). Even at 300 °C, the 
characteristic time associated with cracking reactions is 2300 s - 
far longer than the gas residence time within the reactor and sam¬ 
pling systems. 

The overall mass balance during Prins’ volatile compositions 
experiments ranged from 95% to 98%. However, the total volatile 
and solid yields in Table 1 sum to 100%, therefore it was assumed 
that any volatiles which leaked (<5%wt ini ti a i) had the same average 
composition as those directly measured. 


Table 1 

Willow torrefaction volatile composition experimental conditions and yields. Experimental data from Prins (2005). Modeled yields determined using kinetic parameters from 
Prins (2005) shown in Eqs. (6)-(9). 


Experimentally measured yields 


Modeled yields of pseudo components 


Heating rate 
Heatup period 
Isothermal period 


Water 
Formic acid 
Methanol 
Lactic acid 

Hydroxy acetone 
Carbon dioxide 
Carbon monoxide 
Total volatile yield 
Solid Yield 


VI 


(n n) 

(%Wtj„jtial) 

(%wt ini aal) 

(%Wt iniH ai) 

(%Wtj„jtial) 

(%Wt in iaal) 

(%Wt ini aal) 

(%Wt ini aal) 

(%Wt ini aal) 

(%Wt ini aal) 

(%Wt in iaal) 

(%Wt ini aal) 

(%Wt ini aal) 

(%Wtj„jaal) 

(%Wt ini Hal) 
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2.4. Volatile composition model 

Based on the insights gained from the solid mass loss kinetics 
and volatile composition experiments, a simplified volatile compo¬ 
sition model was developed. Torrefaction has been characterized 
by the fast (lasting 10-60 min) decomposition of the hemicellulose 
and partial decomposition of the cellulose during the first stage fol¬ 
lowed by the gradual (lasting 60-180 min) degradation of the 
remaining cellulose in the second stage. The individual composi¬ 
tions of the pseudo components 1/1 and V2 were each modeled 
with a fixed - but unique - mixture of nine chemical components. 
Since 1/1 and 1/2 are produced in varying proportions depending on 
the temperature and reaction conditions, both the total mass and 
average chemical composition of the volatiles varies with temper¬ 
ature and time. 

The production rate of the nine chemical components (water, 
acetic acid, etc.) depends on the production rates of V 1 and 1/2 
and eighteen model parameters termed compositional coefficients 
{Y a ,\n, VWi, Ycvi -■., etc.) as follows: 

d ™ t * = rviYxyi + r V2 Y xy2 ( 11 ) 

Where d ™' represents the volatile matter production rate in 

(kgxs ') of chemical specie (x = a, b, c . i) and r vl and r V2 are 

the mass production rate of 1/1 and V2 defined by Eqs. (4) and (5), 
respectively. Y xyN is the mass fraction of pseudocomponent VN 

(1/1, V2) composed of chemical component (x = a, b, c .i), in (kg x 

/kg VN). The nine subscripts (a, b, c,.. ., i) refer to acetic acid, water, 
formic acid, methanol, lactic acid, furfural, hydroxyacetone, carbon 
dioxide and carbon monoxide, respectively. 

The chemical compositions of the pseudo components VI and 
V2 might be expected to be dependent on temperature. However, 
given the limited temperature range of the above volatile compo¬ 
sition experiments between 230 and 300 °C, the individual compo¬ 
sitions of 1/1 and 1/2 were assumed constant within this range. 
Moreover, stage one volatiles (1/1) are expected to represent pri¬ 
marily (but not solely) hemicellulose decomposition products 
and stage two volatiles (1/2) should be representative of cellulose 
decomposition products. 

To summarize, the composition of V 1 and 1/2 are modeled with 
unique chemical mixtures of nine species which do not vary with 
temperature (between 230 and 300 °C) and are fitted to experi¬ 
mental data. 

2.4.1. Volatile composition parameter fitting 

The model required fitting 18 parameters. These parameters - 
referred to as compositional coefficients- represented the mass 
fraction contribution to 1/1 and 1/2 of the nine volatiles (acetic acid, 
water, formic acid, etc.). Fitting these parameters required three 
steps: 

(1) The estimation of the amounts of pseudocomponents VI and 
1/2 which would be expected during the five volatile compo¬ 
sition experiments (see Table 1 ) through kinetic simulations. 

(2) The formulation of an inverse problem which related the 
experimentally measured total volatile yields to the kineti- 
cally simulated yields of VI and 1/2. 

(3) The solution for optimal set of 18 compositional coefficients 
through least squares regression. 


2.4.1.7. Kinetics simulation. The evolution of pseudocomponents is 
obtained by integrating Eqs. (1)—(5). By providing appropriate ini¬ 
tial conditions and imposing a particle temperature profile repre¬ 
sentative of the experimental conditions, the cumulative mass 
yields of A, B, C, VI, V2 were obtained as a function of time. The 


simulated VI and V2 yields under torrefaction conditions which 
parallel those in experiments #1-5 are shown in Table 1. The total 
volatile yields Y vtot and solid mass yield Y s are evaluated from the 
simulated yields of VI and V2: 

Yvtot = Y v , + Y V2 (12) 

y s = i-Vet,* = y A + y B + yc (13) 

where Y x represents the yield of X (A B, C, VI, V2, Vtot, S ) in units of 
(kg X/kg initial biomass). 

The kinetics simulation results indicate that pseudocomponent 
VI comprises the vast majority of the total volatile yield under all 
five experimental conditions. V2 comprises between 3.5% and 25% 
of the total volatile yield while VI comprises between 75% and 
96.5%. Because V2 is associated with volatile release during the 
slower second stage, higher torrefaction temperatures and longer 
residence time lead to an increased ratio of V2 relative to VI. 

2.4.1.2. Formulation of the inverse problem. The compositional coef¬ 
ficients (y a ,vi, YW, Y cVl .. etc.) nine describing VI and nine 

describing V2 are used to relate Yvi and Yvi (obtained from the ki¬ 
netic simulations) to those experimentally measured yields of the 
nine chemicals. This relation between the kinetically simulated VI 
and V2 yields and the nine experimentally measured total chemical 
volatile yields is summarized by the following system of equations: 

Y n vl *Y m + Y n V2 *Y xy2 = Y n x (14) 

where Y" is the experimentally measured total yield of volatile spe¬ 
cies (x = a, b, c,..., i) during experiment number (n = 1, 2, 3,.., 5) ex¬ 
pressed as a mass fraction of the initial solid mass (kgx/kg initial 
biomass) (see Table 1). Y% is the simulated yield of VN (VI, V2) 
for experiment (n=l, 2, 3,...,5) expressed as a mass fraction of the 
initial solid mass (Kg VN/kg initial biomass) (see Table 1). This is 
represented in matrix form by the system: 


Yjq y' V2 

y vi Y v2 
Yvi y3 v2 

[Y a ,vi Yfc.vi ... Y m y 

Y\ Y\ ... Yj 
Y 2 a Y\ ... Y] 

Y 3 a Y\ ... Y f 

Yvi Y 4 V2 
.Yvi Y 5 V2 

[Y a y2 Y biV2 ... Y iy2 J 

Yl Y 5 b ... Yf 


where the first matrix on the left-hand side has dimensions of 5 
rows by 2 columns and is composed of the VI and V2 yields gener¬ 
ated by the kinetics simulation, the second matrix has dimensions 
of 2 rows by 9 columns and comprises of the 18 compositional coef¬ 
ficients, and the third matrix has dimensions of 5 rows by 9 
columns and consists of the experimentally reported yields of the 
nine chemical species. 

2.4.1.3. Least squares solution. The 18 compositional coefficients 
were solved iteratively until the least squares solution was found. 
Additional constraints include: 


Ev = i 

(15) 

= i 

(16) 

Y x yi,Y x ,V2 > o 

(17) 


These constrained the least squares solution of 18 composi¬ 
tional coefficients to be physically meaningful by ensuring non¬ 
negative values. Also, since the compositional coefficients are mass 
fractions, they are constrained to sum to 1. 








464 


R.B. Bates, A.F. Ghoniem/Bioresource Technology 124 (2012) 460-469 


2.5. Solid composition model 

The composition of A is fixed and known from the ultimate 
analysis of the raw biomass. The compositions of VI and V2 are 
fixed and defined by the previously described volatile composition 
model. The evolution of the average solid product composition de¬ 
pends only on: 

- The initial solid product composition (i.e. ultimate analysis of 

raw biomass) 

- The amount and composition of the volatiles which have left 

the solid product (i.e. amount and composition of VI, V2). 

However, reaction thermochemistry depends on the composi¬ 
tion of products and reactants generated at the instant at which 
the reactions occur- not on the average composition of products 
and reactants which exist in the system. This requires the compo¬ 
sition of the solid individual pseudo components (B and C) which 
are reacting. 

In summary, two types of solid product composition were 
relevant: 

1) The average solid product composition described in 
Section 2.5.1 

2) The individual compositions of B, C which are reacting/form¬ 
ing at an instant in time. The kinetic analysis regarding the 
instantaneous product distribution is described in Sec¬ 
tion 2.5.2. Section 2.5.3 relates these fractional yields to 
the composition of B and C. 

2.5.1. Average solid product composition 

Combining mass conservation with the kinetics rate expres¬ 
sions, the rate of change of the carbon, hydrogen, oxygen, nitrogen, 
and ash content (in kg) of the solid product is written in terms of 
the composition and formation rates of VI and V2: 

= -r v ,Y jyl - r V2 Yjy 2 (18) 

where, m s is the mass of the solid (in kg). 

Y jf x denotes the mass fraction of pseudo component X (A B, C, VI, 
V2) composed of element (j = C, H, 0, N, Ash) on a dry basis and r vl 
and Ty 2 are previously defined reaction rates (Eqs. 4 and 5, 
respectively). 

2.5.2. Instantaneous fractional yields 

In the present mass loss kinetics mechanism, the pseudocompo¬ 
nents are produced in two stages each with parallel reactions. In 
each stage, decomposition occurs first-order with respect to the 
reactant (A and B for stage 1 and 2, respectively). The dependence 
of the product distribution on temperature is defined by the 
instantaneous fractional yields. Based on this product distribution 
analysis, the composition of B and C can then be defined from mass 
balance. 

The instantaneous fractional yields are positive, dimensionless, 
and defined by the formation rate of product divided by the 


decomposition rate of the reactant: 


„ r B ,i _ kj*m A _ ki 

P ~ -r A ,i ~~ (ki + k vl )* m A ki + k vi 

(19) 

r v 1,1 kvi * m A ki 

-r A ,\ (ki + ky t )*m A k, i k vl 

(20) 

r c , 2 k 2 *m B k 2 

-r B ,2 (k 2 + k V2 ) *m B k 2 + k V2 

(21) 


Ty 2,2 _ ky2 *m B _ k V2 

-r B .2 ~ (k 2 + k n )*m B ~ k 2 + k V2 


(22) 


where, rx, N represents the rate of change of mass (kgjXs ’) of 
pseudocomponent X, (A, B, C, VI, V2) in step N (1, 2) k N are the 
Arrhenius rates with units of 1 /s. Using the instantaneous fractional 
yields we can express the two step kinetic mechanism as: 


A 


[}B + vVl 


(23) 


B^yC+£,V 2 (24) 

where /?, v, y, f are the instantaneous fractional yields. These are 
interpreted in the following way: For each gram of A which is react¬ 
ing at temperature T, /?, grams of B and v grams of VI are formed. 
Likewise, y and £ represent the mass fraction of B transforming into 
C and V2, respectively. By definition, the sum of p and v and the sum 
of equals y and c equals 1. 


2.5.3. Composition of B and C 

From the instantaneous fractional yields defined previously, it is 
possible to define the composition (ultimate analysis) of B and C. 
The composition of B and C as a function of the instantaneous frac¬ 
tional yields and the compositions of A, VI, and V2: 

Y iB = (Y jA -V*Y m )/fl (25) 


y j .c = (Y j , B -k*y i .V2)/y (26) 

where, [S, v, y, { are the instantaneous fractional yields, and Yj tX 
denotes the mass fraction of pseudo component X (A, B, C, VI, V2) 
composed of element (;'=C, H, O, N, Ash) on a dry basis. 


3. Results and discussion 

3.1. Fitted compositional coefficients 

The 18 compositional coefficients obtained by fitting are 
summarized in Table 2. 

The results show that VI is composed primarily of acetic acid, 
water, and carbon dioxide. V2 is composed entirely of condensable 
volatiles including lactic acid, methanol, and acetic acid, water, 
hydroxyacetone, and formic acid. 

These results show agreement with several experimentally ob¬ 
served trends in the volatiles produced during wood pyrolysis. For 
example, a zonal description of the products of wood pyrolysis was 
offered by Gullu and Demirba? (2001). Zone A occurs at tempera¬ 
tures below 200 °C and results in water, carbon dioxide, formic 
acid, acetic acid, glyoxal. Zone B occurs at temperatures between 
200 and 260 °C and results in water, carbon dioxide, formic acid, 
glyoxal, and carbon monoxide. Zone C occurs at temperatures be¬ 
tween 260 and 500 °C and results in methane, formaldehyde, for¬ 
mic acid, acetic acid, methanol, hydrogen. In Zone D, which 
occurs above 500 °C, the final char product is formed. The fitted 
composition of VI agree with that of zone A and B (200-260 °C), 
while the composition of V2 agree with that of Zone C (260- 
500 °C). 

The fitted results also agree with pyrolysis of pure lignocellu- 
lose components. According to Demirba? (2007) the condensable 
products of xylan (hemicellulose) degradation consists primarily 
of eight products including water, methanol, formic acid, acetic 
acid, propionic acids, hydroxyacetone (C 3 H 6 02), and l-hydroxy-2- 
butanone (C 4 H 8 0 2 ). Volatile products of the slow pyrolysis 
(5 °C/min) of galactoglucomannan- the primary type of hemicellu¬ 
lose contained in softwoods- mainly includes water, carbon diox¬ 
ide, carbon monoxide and acetic acid (Aho et al„ 2008). Fast 
pyroylsis ash free hemicellulose derived from switchgrass at 
500 °C result in similar products including carbon dioxide 
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Table 2 

Fitted VI and V2 composition (% mass) through least squares minimization. 
Calculated carbon, hydrogen, and oxygen content. Nitrogen and ash content assumed 
to be negligible in volatiles. 


Chemical Component 


Percentage 

VI V2 


Water 
Formic acid 
Methanol 
Lactic acid 
Furfural 

Hydroxy acetone 
Carbon dioxide 
Carbon monoxide 


Hydrogen 

Oxygen 


y ttV1 = 14.8% 

=48.i% 
y c vi = 5.3% 

y d , vi = 4.2% 

y e , « = 1.3% 

Y fv , = 1.1% 

y s vi - 0.6% 

y d vi = 20.4% 
y t vi = 4.2% 
100% 

y C vi = 18% 
y H ,vi = 7% 
y 0 ,vi = 75% 
y N !v.=o% 
w,=o% 
100% 


y dV 2 = i6.i% 
y 6 ,V2 = 7.6% 
y tV 2=5.i% 

y d v2=30.1% 

y e ,v2 = 31.3% 
y/.v2 = 0.0% 
y s v2 = 9.7% 
y d v2 - 0.0% 
y,-,v2 = 0.1% 
100% 

y C v2 = 36% 
Y H ,V2 = 9% 
Yoy2 - 55% 
y«!vt=o% 

Wi-0% 

100% 


(18.8 wt%), water (15.1 wt%), formic acid (8.4 wt%), and char 
(11 wt%), hydroxyacetone (3 wt%) and other dehydration products 
including dianhydro xyloses(8 wt%) (Patwardhan et al., 2011). 
Lower temperatures (300 °C) and the addition of minerals (alkali 
and alkaline salts) caused increased char yields. The mechanism 
of acetic acid formation from hemicellulose is the elimination of 
acetyl groups from the xylose (C 5 H 10 05) unit (Demirba§, 2000). 
As expected, VI is highly representative of these products. 

Acid catalyzed hydrolysis of hemicellulose results in a mixture 
of pentose (C5) and hexose (C6) sugars and the subsequent dehy¬ 
dration of the pentose sugars results in furfural (C 5 H 4 0 2 ) (Demi- 
rba§, 2000). Therefore, the presence of furfural in VI (1 wt%) and 
absence of furfural in V2 provides further evidence for the proposal 
that VI represents hemicellulose decomposition products. 

The fast pyrolysis of cellulose primarily results in levoglucosan 
(C 6 Hi 0 O 5 ) and other anhydrosugar derivatives (Lu et al., 2011). 
However, levoglucosan is unstable and experiences subsequent 
decomposition and reactions to form volatiles including water, 
methanol, acetic acid, acetone, phenols, (Demirba§, 2000). Formic 
acid, acetic acid, and hydroxyacetone, glycoaldehyde, and phenols 
have been attributed to the ash catalyzed decomposition of levo¬ 
glucosan (Ponder et al., 1992; Patwardhan et al., 2010). Therefore 
the composition of V2 is representative of the pyrolysis of cellulose 
and its anhydrosugars. 

The fitted compositions of pseudo components VI and V2 sup¬ 
port the original assumption that their compositions are represen¬ 
tative of hemicellulose (zone A, B) and cellulose (zone C) pyrolysis 
products, respectively. This strong coupling of the evolution of bio¬ 
mass pyrolysis products to constituent lignocellulosic components 
agrees with the literature (Demirba§, 2007; Aho et al., 2008). 


3.2. Effect of temperature on product distribution 

The instantaneous fractional yields calculated based on the wil¬ 
low kinetics parameters are shown in Fig. la as a function of tem¬ 
perature. Each gram of B that reacts is converted into 0.59 g V2 and 
0.41 g of C regardless of temperature, y and \ are constant due to 
the equal activation energies of k 2 , k v i of that step. However, fi 
and v are temperature dependent due to the unequal activation 
energies of k^, k m . Although the vast majority of A converts into 
the intermediate product (B), increasing amounts of VI (relative 
to B) are produced at higher temperatures. For example, at 
220 °C, each gram of A which is reacting forms 0.9 g of B and 


(a) 

2 1 



200 210 220 230 240 250 260 270 280 290 300 
Temperature (°C) 



0.1 

0.0 


200 210 220 230 240 250 260 270 280 290 300 
Temperature (°C) 

(C) 1 

o 08 
0.6 
§■ 0.4 
>-" 0.2 
0 

200 210 220 230 240 250 260 270 280 290 300 
Temperature (°C) 

Fig. 1. (a) Willow torrefaction instantaneous fractional yields (-) versus final 
temperature (°C). Composition Yp of pseudocomponents B and C shown in (b and c), 
respectively in units of (kgj/kg,- dry basis j = C, H, O i = B, C) versus final temperature 
(°C) as defined by Eqs. 25 and 26. 



0.1 g of VI. However at 300 °C, A degrades to form only 0.7 g of B 
and 0.3 g of VI. This demonstrates the important effect of temper¬ 
ature on the product distribution: higher temperatures shift the 
product distribution towards greater volatile product formation. 

The competitive nature of the first step of this kinetic mecha¬ 
nism is representative of the pyrolysis characteristics of pure xylan 
(hemicellulose) demonstrated during TGA mass loss kinetics 
experiments by Di Blasi and Lanzetta (1997) and van der Stelt 
(2011). In these experiments, pure xylan was subjected to a 10- 
70 °C/min heating rate until reaching a prescribed maximum tem¬ 
perature (between 220 and 300 °C) after which it was held at this 
temperature for up to 120 min. The maximum mass loss (i.e. vola¬ 
tile formation) during this xylan pyrolysis exhibited temperature 
dependence. For example, the final volatile yield of ranged from 
42 to 56%wt initial after 120 min depending on the final temperature 
(240-300 °C, respectively) (van der Stelt, 2011 ). In other words, 
higher temperatures favors the formation of volatiles (relative to 
solid product) during hemicellulose pyrolysis. 

Cellulose pyrolysis also demonstrates a competitive behavior 
between volatile and char formation. In TGA experiments of cellu¬ 
lose between 240 and 300 °C in (Cho et al., 2010), the final volatile 
yield ranged between 70 and 90%wt init iai (increasing linearly with 
temperature). 
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(b) 



Time (minutes) 


Fig. 2. Modeled solid yield composition (a) and volatile yie 

This suggests that the mechanism of volatile formation during 
the first step (A -» B + VI ) is primarily (but not solely) attributed 
to hemicellulose decomposition, with an increasing contribution 
from cellulose at higher temperatures. Based on the results of 
TGA mass loss ldnetic studies of cellulose, hemicellulose and lignin 
torrefaction for one hour at (230, 260, and 290 °C) (Chen and Kuo, 
2011) also drew this conclusion. Hemicellulose and cellulose 
decomposition were classified as being as strongly active for above 
temperatures of 260 and 290 °C, respectively. In summary, not only 
is the rate of torrefaction strongly more rapid at higher tempera¬ 
tures, but the product distribution shifts towards increased volatile 
formation. 


3.3. Effect of temperature on solid product composition 

Substituting the instantaneous fractional yields, the initial ulti¬ 
mate analysis of willow (Y jA = 0.472, 0.061, 0.451, 0.003, 0.013, 
j = C, H, O, N, Ash) (Prins, 2005) and the composition of VI, V2 
(Yjyi and Yj iV2 , (see Table 2) into Eqs. 25 and 26 results in Fig. lb 
and c which summarize the instantaneous composition of B and 
C as a function of temperature, respectively. At 300 °C (the upper 
temperature limit for torrefaction), C, which represents char, is 
composed almost entirely of carbon (94%). At the 200 °C (the lower 


[ and composition (b) (% weight of initial) versus time (min). 

temperature limit for torrefaction) B has a very similar composi¬ 
tion to A (non torrefied willow). In both B and C, as temperature in¬ 
creases the mass fraction of carbon increases while those of 
hydrogen and oxygen decreases. This trend of reduced hydrogen 
and oxygen mass fraction in the char with higher peak tempera¬ 
ture has been noted in Neves et al. (2011) which performed a 
structured collection of pyrolysis product composition data over 
a range of pyrolysis peak temperatures (200-1000 °C). 

Based on the decomposition of the pure components, it can be 
concluded that B chemically represents an intermediate solid con¬ 
taining the charred remnants of hemicellulose decomposition 
along with partially reacted cellulose and weakly degraded lignin 
fraction. C represents a final charred product containing com¬ 
pletely pyrolyzed hemicellulose and cellulose with a partially 
decomposed lignin fraction. Thus, the application of higher torre¬ 
faction temperature results in a more carbon dense solid product 
than lower torrefaction temperatures. 

3.4. Effect of residence time on solid and volatile product composition 

The effect of residence time on the average solid product com¬ 
position is obtained through integration of Eqs. 11 and 18. For this 
integration, the solid product composition of willow (Prins, 2005) 
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Mass loss (kg/kg initial) 


Fig. 4. Yield of gaseous and liquid (condensable) volatiles during willow torrefaction versus mass loss. Experimental data points from (Zanzi et al„ 2002). 


was used and the initial solid temperature was 200 °C ramping at 
10 (°C/min) to the final temperature of 300 °C, which is represen¬ 
tative of TGA conditions. 

Fig. 2 illustrates both the overall (solid mass loss) as well as the 
volatile product composition with time. This composition profile is 
in line with the expected de-oxygenation of the solid product that 
occurs due to the decarboxylation and dehydration reactions (evi¬ 
denced by the high carbon dioxide and water content present in 
the volatile products). Also of importance is the relatively high car¬ 
bon yield in the solid product throughout the process. After 15 min 
of residence time, the torrefied product loses 40% of its original 
oxygen content but only 11% of the carbon content through the 
evolution of volatiles. At 40 min this increases to 69% and 26%, 
respectively. 

Initially, the volatile composition is composed mostly of VI- 
characterized by high water, acetic acid and carbon dioxide con¬ 


tent). As the torrefaction proceeds, V2 is produced in greater quan¬ 
tities resulting in increasing yields of lactic acid, methanol, and 
acetic acid. Though the fitted compositions are fixed, the yield 
and composition of total volatiles varies because the pseudocom¬ 
ponents (VI and V2) are produced in varying amounts and propor¬ 
tions as residence time increases. As the second stage begins to 
predominate, the carbon content of the volatiles increases. More¬ 
over, the rate of mass loss also becomes slower. As a result of these 
two effects, the overall rate of de-oxygenation during the second 
step is much slower than the first step. The behavior suggests that 
limited de-oxygenation is gained by increasing residence time 
torrefying biomass beyond the time required to complete the first 
stage. Similar conclusions regarding the optimal residence time 
were drawn by Prins (2005). 

The shift in composition of volatiles is demonstrated by the 
time derivative of the volatile yields shown in Fig. 3. The present 
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Mass loss (kg/kg initial) 


Fig. 5. Yield of element j retained in solid product per unit mass of element j in initial biomass, Y SJ /Yaj j = C, H, 0 on a dry ash free basis. Data points frc 
et al. (2008) and Medic et al. (2011). 


i Prins (2005), Arias 


results indicate a peak rate of volatile release at approximately 
10 min which is dominated by water, acetic acid, and carbon diox¬ 
ide production. From 20 min onwards, the volatile release is pri¬ 
marily composed of lactic acid and methanol. These results share 
qualitatively similar trends with the FTIR spectra of reed canary 
grass at 290 °C published in Bridgeman et al. (2008) which also 
demonstrates a single peak in volatile production occurring as 
the sample reaches the final temperature. Water and carbon diox¬ 
ide were found to be the primary volatile constituents during this 
initial peak in release, agreeing with the present model results. 

3.5. Validation with experimental data 

Model predictions for the cumulative yield of non-condensable 
and condensable volatiles versus solid mass loss are compared 
with experimental data on willow torrefaction reported by Zanzi 
et al. (2002) in Fig. 4. The model predicts that the liquid (condens¬ 
able components such as water, acetic acid, methanol, etc.) 
volatiles comprise 70-80% of the total volatile yield. This agrees 
with the data from Zanzi which observed that 80-90% of the 
volatiles were condensable. Moreover, the C0 2 to CO ratio 
predicted by the model (4.85) is within experimentally measured 
range between 3.3 and 7.7 (Zanzi et al., 2002). Thus the model 
satisfactorily predicts the distribution of condensable and non¬ 
condensable volatile products. 

Using results shown in Fig. 2 we can evaluate the fraction of 
carbon, hydrogen, and oxygen retained by the solid product versus 
the mass loss. In Fig. 5, model predictions for these quantities are 
compared with experimental data from Prins (2005), Arias et al. 
(2008) and Medic et al. (2011) who measured the initial and final 
ultimate analyses of torrefied willow, com stover, and eucalyptus 
samples, respectively over a range of torrefaction temperatures 
between 200 and 300 °C and residence times up to 45 min. 

The model predictions agree with the experimentally measured 
data which demonstrate the carbonizing effect of torrefaction on 
the residual solid product: the retained fraction of carbon is always 
higher than the retained fractions of oxygen and hydrogen. These 
results suggest that the present model satisfactorily characterizes 
the evolution in solid product composition for a variety torrefac¬ 
tion conditions and feedstocks other than willow. However, the 
model predictions appear to over-predict the solid product carbon 


yield and under predict the oxygen yield at higher levels of mass 
loss. 

One reason why the model may overpredict the carbon yield is 
the assumption regarding the unaccounted or leaked volatiles dur¬ 
ing the volatile composition experiments. It was assumed that any 
leaked volatiles had the same composition as those measured; 
however, certain carbon rich volatile products could have com¬ 
prised a greater portion of the unnaccounted volatiles than those 
measured in the HPLC and GC. The condensation of compounds 
within the tubing and fittings might cause this. If such experimen¬ 
tal effects had occurred, the model would tend to underpredict the 
carbon content of the volatile products and therefore over-predict 
the solid product carbon yield. 

Another reason is that the kinetics and composition parameters 
do not explicitly account for variations and interactions of the lig- 
nocellulose components, mineral content, and heating rate effects 
which can have important effects on volatile yields Patwardhan 
et al. (2010) and Wang et al. (2011). The experimental data repre¬ 
sented in Fig. 5 cover a range of feedstocks: corn stover has a dras¬ 
tically different lignocellulose composition and ash content from 
willow and would certainly have unique solid mass loss kinetic 
parameters. The composition of volatiles could be expected to be 
different as well. 

Despite feedstock specificity of the model parameters, the 
methodology for fitting the parameters is feedstock general: 
similar mass loss kinetics (using TGA) and volatile composition 
experimental data (using GC/HPLC) should be applied to other 
biomass feedstocks. 


4. Conclusion 

A model of the evolution of volatile and solid product composi¬ 
tion was proposed and the parameters for the model have been fit¬ 
ted from published solid mass loss kinetic and volatile composition 
experiments for kinetically limited willow torrefaction experi¬ 
ments between 230 and 300 °C. 

The model describes the rate of release and composition of vol¬ 
atiles in terms of nine identifiable components enabling the com¬ 
position of the solid product to be predicted through elemental 
balances. The model predictions of the solid product residue com- 
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position and volatile release agree with published experimental re¬ 
sults which were not used to derive the parameters. 
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